Triplet contribution to the Josephson current in the nonequilibrium 
superconductor / ferromagnet / superconductor junction 
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The Josephson current through a long s-wave superconductor/weak ferromagnet/s-wave super- 
conductor weak link is studied theoretically in the regime of nonequilibrium spin-dependent occu- 
pation of electron states in the ferromagnetic intelayer. While under the considered nonequilibrium 
conditions the standard supercurrent, carried by the singlet part of current-carrying density of 
states, practically is not modified, the additional supercurrent flowing via the triplet part of the 
current-carrying density of states appears. Depending on voltage, controlling the particular form of 
spin-dependent nonequilibrium in the interlayer, this additional current can enhance or reduce the 
usual current of the singlet component and also switch the junction between 0- and 7r-states. 

PACS numbers: 74.45. +c, 74.50.+r 



New states have been predicted and observed in 
Josephson weak links in the last years. One of them 
has common origin with the famous LOFF-statei 2 . This 
mesoscopic LOFF-state, which is induced in supercon- 
ductor/weak ferromagnet /superconductor (SFS) Joseph- 
son junction, was predicted theoretically^ and observed 
experimentally 5-5 . In this state Cooper pair acquires the 
total momentum 2Q or — 2Q inside the ferromagnet as 
a response to the energy difference between the two spin 
directions. Here Q oc h/vt, where h is an exchange en- 
ergy and Vf is the Fermi velocity. Combination of the two 
possibilities results in the spatial oscillations of the con- 
densate wave function ^(x) in the ferromagnet along the 
direction normal to the SF interface. ^ s (x) oc cos(2Qa;) 
for the singlet Cooper pair—. The same picture is also 
valid in the diffusive limit. The only thing we need to 
add is an extra decay of the condensate wave function 
due to scattering. In the regime /i> |A|, where A is a 
superconducting order parameter (OP) in the leads, the 
decay le ngth is equal to the magnetic coherence length 
£f = y D/h, while the oscillation period is given by 
2tt^f- Here D is the diffusion constant in the ferromag- 
net, H = 1 throughout the paper. 

The presence of an exchange field also leads to the for- 
mation of the triplet component of the condensate wave 
function in the interlayer. In the case of a homogeneous 
exchange field only the component with zero spin pro- 
jection on the field direction S z = is induced. Com- 
bining the two pairs with the total momenta 2Q and 
— 2Q into the triplet combination we see that the in- 
plane (S z — 0) triplet condensate wave function com- 
ponent ^t(x) oc sin(2Qx), that is oscillates in space with 
the same period as the singlet one, but is shifted by ir/2 
with respect to it. Here we do not discuss the other triplet 
components with S z = ±1, which are typically induced 
in case of inhomogeneous magnetization 10 . 

The energy spectrum of the superconducting correla- 
tions is expressed in a so-called supercurrent- carrying 
density of states (SCDOS)^ - — . This quantity represents 
the density of states weighted by a factor proportional to 
the current that each state carries in a certain direction. 



Under equilibrium conditions the supercurrent can be ex- 
pressed via the SCDOS as^ 



J deN J (e)t&nhe/2T 



(1) 



where e stands for the quasiparticle energy, <p(s) = 
tanhe/2T is the equilibrium distribution function and 
Nj(e) is SCDOS. In the presence of spin effects SCDOS 
becomes a matrix 2 x 2 in spin space and can be repre- 
sented as Nj — Nj. s + Njjtr, where Ui are Pauli matrices 
in spin space. The scalar in spin space part of SCDOS 
Nj i8 is referred to as the singlet part of SCDOS in the 
paper and the vector part Njj is referred to as the triplet 
part. Nj.t is directly proportional to the triplet part of 
the condensate wave function. It is well known that the 
spin supercurrent can not flow through the singlet super- 
conducting leads. Therefore, Njj does not contribute to 
the supercurrent in equilibrium. Having in mind that the 
triplet part of SCDOS is even function of quasiparticle 
energy, one can directly see that this is indeed the case. 

Under nonequilibrium conditions one can change the 
value of the critical Josephson current through the 
junction and even realize the 7r-state by manipulating 
the quasiparticle distribution in an interlayer region. 
This effect was predicted theoreticall y 1 ^ 12 and observed 
experimentally 15-17 for a diffusive SNS junction. The 
point is that positive and negative parts of SCDOS give 
energy-dependent contributions to the supercurrent in 
the positive and negative direction. The size and di- 
rection of the total supercurrent depends therefore on 
the occupied fraction of these states, which is analogous 
to the occupation of the discrete Andreev bound states 
in a ballistic system. That is one can obtain negative 
Josephson current response to small phase differences 
and, hence, switch the system into 7r-state by creating 
an appropriate nonequilibrium quasiparticle distribution 
in the weak link region. The combination of the exchange 
field h <C |A| and the spin independent nonequilibrium 
distribution function has been considered as welU^. Un- 
der these conditions the influence of the nonequilibrium 
distribution function is also consists in the redistribution 
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of the quasiparticles between the energy levels. How- 
ever, it was shown that in this limit of small exchange 
fields the combined effect of the exhange field and the 
nonequilibrium distribution function is nontrivial. For 
instance, part of the field-suppressed supercurrent can 
be recovered by adjusting a voltage between additional 
electrodes, which controls the ditribution function. 

In the present paper we investigate the effect of 
nonequilibrium occupation of the supercurrent-carrying 
states on the Josephson current in SFS junction in the 
parameter range |A| <C h <C £/, where e/ is the Fermi en- 
ergy of the ferromagnet. This regime is relevant to weak 
ferromagnetic alloys, which were used for the experimen- 
tal realization of magnetic 7r-junctions. It is shown that 
if the distribution function becomes nonequilibrium and 
spin-dependent, the supercurrent carried by the SCDOS 
triplet component Nj t in the ferromagnet is non-zero. 
The magnitude of this current cotribution j t can be of 
the same order or even larger than the current contri- 
bution j s carried by the singlet component Nj iS . Due 
to the fact that the singlet and triplet components of 
the anomalous Green's function have the same oscillation 
period but shifted in phase by tt/2, jt can increase the 
usual supercurrent, carried by the singlet component of 
SCDOS, or weaken it, or even reverse the sign of the total 
supercurrent, thus switching between and 7r-state. Ex- 
perimentally the most probable way to realize the spin- 
dependent nonequilibrium in the interlayer is to apply a 
voltage to (or to pass the dissipative current through) a 
spin-active material. Then under appropriate conditions 
even quite small voltages should be enough to switch the 
system from to 7r-state and vice versa. 

That is we consider another mechanism of supercurrent 
manipulation by creating a nonequilibrium quasiparticle 
distribution in the interlayer, which cannot be reduced to 
the redistribution of the quasiparticles between the en- 
ergy levels. In principle, the both mechanisms can be re- 
alized in a junction simultaneously. However, in this par- 
ticular study we assume /i > A and the Thouless energy 



£Th = D/d 2 <C h, that is the interlayer length d >• £p. 
As it is shown below, in this regime (and for the case 
of low-transparency SF interfaces) a spin- independent 
nonequilibrium distribution of quasiparticles in the ferro- 
magnet practically does not affect the Josephson current, 
that is the described above mechanism of the critical cur- 
rent reversal by the spin-independent redistribution of 
supercurrent-carrying states population is irrelevant in 
this case. For high-transparency SF interfaces the both 
mechanisms contribute to supercurrent. 

For a quantitative analysis we use the formalism of 
quasiclassical Green-Keldysh functions in the diffusive 
limits. The fundamental quantity for diffusive transport 
is the momentum average of the quasiclassical Green's 
function g{x,e) — (g(pf,x,e)) Pf . It is a 8 x 8 matrix 
form in the product space of Keldysh, particle-hole and 
spin variables. Here x is the coordinate measured along 
the normal to the junction. 

The electric current should be calculated via Keldysh 
part of the quasiclassical Green's function. For the plane 
diffusive junction the corresponding expression reads as 
follows 



+ OC 



de 
8^ 



rTr 4 



T() 



— (g(x,e)d x g(x,e)) K 



(2) 



where e is the electron charge and Rf stands for the re- 
sistance of the ferromagnetic region. (g(x, e)d x g{x, e)) 
is 4 x 4 Keldysh part of the corresponding combination of 
full Green's functions. Tj are Pauli matrices in particle- 
hole space. 

It is convenient to express Keldysh part of the full 
Green's function via the retarded and advanced compo- 
nents and the distribution function: g K = g R (p — *pg A - 
Here the argument (x, e) of all the functions is omitted for 
brevity. The distribution function is diagonal in particle- 
hole space: ip — (p(T + T 3 )/2 + a2'f (72 {tq — t 3 )/2. All the 
matrices denoted by " are 2x2 matrices in spin space 
throughout the paper. In terms of the distribution func- 
tion current ^ takes the form 



J = 



-d 
eRp 



de 
8^ 



:Tr 2 



i?d*$ - g H d x <pg A - f H d x £f A + (g H d x g H + f H d x f H )0 - ^(g A d : 



f A d x f A ) 



(3) 



We assume that the direction of the exchange field h is 
spatially homogeneous and choose the quantization axis 
along the field. In this case the distribution function 
and the normal part g R > A of the Green's function are di- 
agonal matrices in spin space. The anomalous Green's 
functions can be represented as f R ' A = f^' A i<J2 and 
Jr,a _ _i (T2 fR, ^ where f R ' A and f R ' A are diagonal 
in spin space. 



The retarded and advanced Green's functions are ob- 
tained by solving the Usadel equations^ supplemented 
with Kupriyanov-Lukichev boundary conditions at SF 
interfaces^. It is worth to note that we can safely apply 
these boundary conditions to the problem of plane diffu- 
sive junction even for high enough dimensionless conduc- 
tance g of the SF interface. This can be done in spite of 
the fact that they are the linear in transparency approx- 
imation of more general Nazarov boundary conditions 21 . 
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The reason is that the effective number of interface chan- 
nels N ~ dy/l is large and a separate channel trans- 
parency T ~ g(l/d v ) is usually considerably less than 
unity. Here d y is the junction width and I is the mean 
free path. 

Further, the condition d £f allows us to find the 
solution analytically even for an arbitrary SF interface 
transparency and low temperature, that is in the param- 
eter region, where the equations cannot be linearized with 
respect to the anomalous Green's function. We start from 
the completely incoherent junction (that is, consider the 
left and right SF interfaces separately) and then calculate 
the corrections up to the first order of the small param- 
eter exp [— c2/£f] to the Green's functions. Within this 
accuracy the anomalous Green's functions in the vicinity 
of left and right SF interfaces (at x = take the 

form 



fda = Km 



smhQf A e- iax ' 2 + AX R > A (x)e 



iax/2 



fR,A _ _ rR.A / 
J da J da vA 



-X) • (4) 

Here a =t, I (or +1(— 1) within Equations) is the electron 
spin projection on the quantization axis, K = +1(— 1) 
corresponds to the retarded (advanced) functions, a = 
— 1) in the vicinity of the left (right) SF interface 
and x is the order parameter phase difference between 
the superconducting leads. The first term represents the 
anomalous Green's function at the ferromagnetic side of 
the isolated SF boundary and does not enter the following 
results. So, we do not give it explicitly. The second 
term is the first order correction, originated from the 
anomalous Green's function extended from the other SF 
interface: 



y 2 



C )(i 



X R > A (x)=K?' A e 
where K R ' A is determined by the equation 



)/6 



(5) 



(1 



i K a)K R ' A (l-K R ' A 



2 _ Rf^F 

' ~ 4Rgd 



[sinh6f' A (l- 



K R,A* 



cosh 9f ' A 4K?> A {1 + K?* ) (6) 



R g stands for the resistance of each SF interface, which 
are supposed to be identical. sinh0f ,j4 and cosh6f ,A 
originate from the Green's functions at the superconduct- 
ing side of SF interfaces. We assume that the parameter 
(RF£,s/Rgd)(<JF/&s) *C 1, where £ s = y/D/A is the su- 
perconducting coherence length in the leads, of and a s 
stand for conductivities of ferromagnetic and supercon- 
ducting materials, respectively. It allows us to neglect 
the suppression of the superconducting order parameter 
in the S leads near the interface and take the Green's 
functions at the superconducting side of the boundaries 
to be equal to their bulk values. In this case 



coshef' 71 = 
sinhe5 A = 



— me 



vW - ( e + nisy 

— Kl\A\ 



(7) 



where 5 is a positive infinitesimal. 

The SCDOS, entering the current ([3]) takes the form 



N ja (e) = (g R d x g R + f R dJ R 
8tt 2 sin x. 



OpRq 



-Im 



Z?(x 



'd x g A - f d x f J 



-a-)sinh6 6 



(8) 



The nonequilibrium distribution function in the inter- 
layer is proposed to be created by applying a voltage 
along the y direction between two additional electrodes 
Ni, and N t , which are attached to the central part of the 
interlayer. It is supposed that the conductances of 2V&F 
and NtF interfaces are spin-dependent and equal to gb a 
and g ta , respectively. The voltage V ta (Vb a ) between the 
superconducting leads and Nt {Nb) electrode can also be 
spin-dependent. It can be realized, for example, by at- 
taching one (or both) of the electrodes Nb or Nt to a 
strong ferromagnet and applying a voltage between the 
other one and the ferromagnet. 

In order to obtain the distribution function up to the 
first order of the parameter exp [— cZ/£f] we solve the ki- 
netic equation for the distribution function, which is de- 
rived from the Keldysh part of the Usadel equation. The 
boundary conditions to the kinetic equation are also ob- 
tained from the Keldysh part of the general Kypriyanov- 
Lukichev boundary conditions. Further it is assumed 
that |eVt,b*. J < |A|. Under this condition the part of the 
current associated with the first three terms in Eq. ([3]) is 
only determined by the first order correction T, R ' A to the 
anomalous Green's function and the distribution function 
<Pa°\ calculated up to the zero order of the parameter 
exp [-d/Zp]-- 



(-i?d x Cp - g R d x p>g A - f R dJf A )= Nj, 



c3 (0) - w (0) 



(9) 

where Nj a (s) is expressed by Eq. ((5J. 

The distribution function tp^ does not depend on x. 
So, it is convenient to calculate it in the middle of the in- 
terlayer, where disregarding the parameter exp [— cZ/£f] 
means disregarding the superconducting proximity ef- 
fect. Then, by considering Nb/F/Nt junction and ap- 
plying Kypriyanov-Lukichev boundary conditions Nb/F 
and F/Nt interfaces, we come to the following expression 
for <p^ (inelastic scattering processes are not taken into 
account): 



tanh 



s-eVtc 



? (0) _ haLLLL 2T ' 9tg{^F + d v9bg) 
VF{gt a + 9ba) + 2d v9ta9ba 

^ tanh e ~2T b " 9b a {o-F + d y g t(T ) 

OF{gt a + gb a ) + 2d y9ta9ba 



(10) 



<Pa is connected to p^ by the symmetry relation 



(^j°j(e) = — </?j lJ |(— e). We focus on the case gta -C OFjdy 
or gba <C OF/dy, when the y-dependence of the distribu- 
tion function tp^ can be disregarded. 
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Substituting Eqs. ©-© into Eq. ^ we find that the 
Josephson current takes the form 



j = 



-d 
2eR F 



where Nj a (e) is expressed by Eq. and ip^ should 
be taken from Eq. (flOl) . Elecrtical current (11) can be 
divided into two parts: 



Js = Js.cSmx = 
jt =jt, c sin x 



J 

-d f 

eR F J 


= Js^ 

de 
8^ 


" Jt 1 
(4 0) 4 






-d ! 

eR F J 


r de 

8^ 


\^ - 


f4 0) )% 


(12 



where tp Q = (ip f + ip\)/2, ip z = (ip t - <Pi)/2, ip a . z = 
(<p t ± <f>i)/2 = t^(-s), Nj jS = (jy,- t + A^)/2 is the 
singlet part of SCDOS and N j>t = (N jt - N oi )/2 is the 
z-component of the SCDOS triplet part (the other com- 
ponents equal to zero for the considered case of homo- 
geneous magnetization). It is seen from Eq. (fT2"|) that 
Njj gives rise to the additional contribution to the spin- 
less electrical current if the quasiparticle distribution is 
spin-dependent . 

While Eqs. ([Tl" j) - (jT2"]l are valid for arbitrary SF inter- 
face transparency, at first we concentrate on the discus- 
sion of the tunnel limit g = Rf£f/ Rgd <C 1, where 
Eq. (J6j) can be easily solved and the integral over energy 
can be calculated analytically. For the analytical analy- 
sis we choose the most simple form for the distribution 
function Eq. (TTU1) by setting g tu -» and = 0. Then 

<^ 0) = tanh[(e - eV bf )/2T], while ipf ] = tanh[e/2T]. As 
it is demonstrated below, the results corresponding to an- 
other set of parameters, which determines the particular 
form of the distribution function satisfying Eq. (|10j) . qual- 
itatively do not differ from that ones represented here. At 
low temperature T -C \eV b ^\ we obtain for the Josephson 
current the following result 



RfCf sinx I 



AeR 2 g d 



-d/SF 



i— , d 7T 

V27TCOS(- h — ) 

hp 4 



V2 



log 



lAI+eVt 



bf 



IAI 



sin( 



(13) 



The first term in Eq. (|13l) represents the part of the 
supercurrent j s carried by the singlet component of SC- 
DOS. Under the conditions T < \eVt,b\ and |eVt,6* J < 
|A| it is not affected by the fact that the distribution 
function is nonequilibrium, as can be seen from Eq. (1131) . 
The reason is that for a long junction and /i> | A| in the 
tunnel limit the singlet part of SCDOS is concentrated 
in the narrow energy intervals around e = ±|A|, as it is 
illustrated in panel (a) of Fig. [1] This is opposed to the 
cases of diffusive SNS junction with h = 0^— and SFS 
junction with h <C A 1 *, where the singlet part of SCDOS 
is finite and exhibits nontrivial energy dependence in the 
subgap energy region |e| < |A|. Under the conditions 



- (sgn[e - eVbf] + sgn[e + eV bt ] + 2sgn[e])/2, 



(°) , 
<Po +<Pb 

which enters the expression for j s [Ea.(|12p]. practically 
does not differ from its equilibrium value for energy inter- 
vals around e — ±|A|. That is, under the considered con- 
ditions the widely discussed in the literature mechanism 
of supercurrent manipulation by nonequilibrium redistri- 
bution of quasiparticles between energy level o 11 ' 12 ' 15 " — 
does not contribute. 

The second term j t is caused by the triplet component 
of SCDOS and vanishes in the equilibrium \%f — 0. As it 
is seen in panel (b) of Fig. [U the triplet part of SCDOS 
Njj is an even function of energy and has finite value 
in the subgap energy region. So, multiplying it by the 
distribution function 
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(o) 



(sgn[e - eVbf] - sgn[e + eV bt ])/2 , (14) 



one obtains current contribution j t . The absolute value of 
this contribution is roughly proportional to V b f for small 
enough values of this parameter and increases sharply 
when approaches |A|. This behavior is a consequence 
of two facts: (i) distribution function (fT4|) is a constant 
within energy interval [— eV b \, eV b \\ and vanishes outside 
it, (ii) the triplet part of SCDOS has a particular shape 
shown in panel (b) of Fig. [TJ 

To calculate the Josephson current for the case of ar- 
bitrary transparency of SF interfaces one needs to solve 
Eq. © numerically and make use of Eq. (fTTj) . The re- 
sulting curves as functions of eV&t are plotted in Fig. [2] 
Panel (a) shows the current for low enough dimcnsionless 
conductance of SF interface g = 0.1, while panel (b) rep- 
resents the case of highly transparent interface g = 3.0. 
Different curves correspond to different lengths d of the 
junction. In dependence on V b f the current can be en- 
hanced or reduced with respect to its value at V b ^ = 0. 
If the length of the equilibrium junction is not far from 
the 0-7T transition, then small enough voltage can switch 
between the states. Separate plots of the current contri- 
butions j s and j t are represented in panel (a) of Fig.[3]for 
the low-transparency junction with g = 0.1 and in panel 
(b) for the high-transparency junction with g = 3.0. It 
is seen, that for low-transparency junction the current 
and its separate contributions j s and j t behave just as 
described by the tunnel limit discussed above. While the 
tunnel limit qualitatively captures the essential physics 
for the high-transparency junction as well, there are some 
new features, which are discussed below. 

The singlet part of SCDOS for the case of high- 
transparency junction is plotted in panel (c) of Fig. [T] 
It is clearly seen that as distinct from the limit of tun- 
nel junction, it is not only concentrated around super- 
conducting gap edges, but is finite in the whole subgap 
region. It results in the sensitivity of j s to the fact that 
the distribution function ip^ +if^ is nonequilibrium. In 
other words, upon increasing of the SF-interface conduc- 
tance the mechanism of the current controlling based on 
the redistribution of the quasiparticles between energy 
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FIG. 1. The singlet and triplet parts of SCDOS in dependence on the quasiparticle energy. The left column corresponds to the 
low-transparency limit g — 0.1, while the right column represents the high-transparency case g = 3.0. 
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FIG. 2. The critical Josephson current as a function of 
e\4.f/A. The different curves correspond to different lengths 
d of the junction, which are measured in units of £f- 



levels starts to contribute. It is seen from panel (b) of 
Fig.[3]that the difference j s ,c(Vb-f) — j s ,c(Hf = 0) consid- 
erably grows when |eV&f| approaches |A|. It is worth to 
note that the absolute value of j SjC is always reduced by 
this mechanism. This fact leads to the reduction of the 
absolute value of the high-transparency total supercur- 
rent (see Fig. HJb)) at eVb-f —> A. This behavior should 
be compared to the tunnel limit, where the dependence 
of js.c on Vb\ is negligible and absolute value of the total 
supercurrent only increases for eVb\ — > A due to growing 
contribution of j^ c . 

Now we discuss the dependence of the obtained results 
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FIG. 3. Separate plots of the current contributions j 3tC 
(dashed line) and jt, c (dotted line) together with the total 
critical current j StC + jt,c (solid line) as a function of eVbf/A. 
Panel (a) represents the low-transparency case and panel (b) 
corresponds to high-transparency case. 



to the particular shape of the distribution function. It 
is illustrated in Fig. |H Four different examples of the 
distribution function z-component (p® + (p^ 1 , satisfying 
Eq. (|10j) are shown in the right panel. Due to symme- 
try relations this combination is always an even function 
of energy. The scalar part of the distribution function 
is not plotted in the figure because it does not influence 
the supercurrent in the tunnel limit. The correspond- 
ing plots of the critical current in the low-transparency 
limit are represented in the left panel. It can be con- 
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FIG. 4. Dependence of the critical current on the partic- 
ular shape of the distribution function. The right panel 
represents four different examples of the distribution func- 
tion vector part. The curves marked by the numbers from 
1 to 4 correspond to the following sets of parameters in 
Eq. HO]): 1. g ta ->• 0, V H = 0; 2. g ttl g H -> 0, eV U = 0.8A; 
3- 9tt,9u ->■ 0, eV ti = -0.8A and 4. g n = 0.1a F /d y , 
gts, = 4a F /d y , g b f = 3a F /d y , g b ,i = 0.2a F /d y , Vb± = 0.514^, 
eVt-f = 0.35A, eVti — 0.7A. The respective critical cur- 
rent plots in dependence on eVbf, calculated for the same 
sets 1-4 of the distribution function parameters are shown in 
the left panel. The dashed vertical line marks the position 
eVbf = 0.2A, for which the distribution functions are calcu- 
lated. 



eluded that, while the quantitative value of the critical 
current depends on the particular choice of the distribu- 
tion function, this choice has no qualitative effect on the 
supercurrent behavior, as it was already pointed before. 

In conclusion, we have studied the effects of nonequilib- 
rium spin-dependent electron distribution in a weakly fer- 
romagnetic interlayer on the Josephson current through 
SFS junction. It is shown that the nonequilibrium spin- 
dependent electron distribution gives rise to the super- 
current carried by the triplet component of SCDOS. De- 
pending on voltage, controlling the particular form of 
spin-dependent nonequilibrium in the interlayer, this ad- 
ditional current can enhance or reduce the usual current 
of the singlet component and also switch the junction 
between 0- and 7r-states. 
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